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An approximate method is suggested to obtain analytical expressions for the eigenvalues and eigen- 
functions of the some quantum optical models. The method is based on the Lie-type transformation 
of the Hamiltonians. In a particular case it is demonstrated that E x e Jahn- Teller Hamiltonian 
can easily be solved within the framework of the suggested approximation. The method presented 
here is conceptually simple and can easily be extended to the other quantum optical models. We 
also show that for a purely imaginary coupling the E x e Hamiltonian becomes non-Hermitian but 
P(Jo-symmetric. Possible generalization of this approach is outlined. 

PACS numbers: 03.65.Fd, 42.50.Ap 



INTRODUCTION 



It is well known that the rotating wave approximation (RWA) is a useful method in determination of the eigenvalues 
and associated eigenfunctions of the various quantum optical Hamiltonians. The approximation gives accurate results 
when the frequency associated with the free evaluation of the system is essentially bigger than the transmission 
frequencies induced by the interaction between subsystem or external source. In quantum physics the application 
of the RWA usually leads to symmetry breaking: the representation space of the whole system is then divided into 
invariant subspaces, which strongly simplifies the mathematical complexity of the problem and usually provides the 
exact solution of the Hamiltonian. 

The simplest model which describes a two-level atom interacting with a single mode cavity field is the Jaynes- 
Cummings (JC) model Jj. A considerable attention has been devoted to the interaction of a radiation field with 
atoms since the paper of Dicke Q . Such system is commonly termed as the Dicke model. In spite of its simplicity, 
the whole spectrum of the Dicke Hamiltonian can not be obtained exactly and usually it has been treated in the 
framework of RWA. Besides its solution with RWA, in some papers an attempt is made to go beyond the RWA j^. 
The continual integration methods are based on variational principles. The perturbative approach leads to more 
complicated mathematical treatments and the theory converges only for a certain relationship between parameters 
of the Hamiltonian. In a more recent study, Klimov and his co-workers have developed a general perturbative 
approach to quantum optical models beyond the RWA, based on the Lie-type transformation. 

The Jahn Teller (JT) interaction [3] is one of the most fascinating phenomena in modern physics and chemistry, 
providing a general approach to understanding the properties of molecules and crystals and their origins. This 
phenomena has inspired of the most important recent scientific discoveries, such as the concept of high temperature 
superconductivity. The JT interaction is an example of electron-phonon coupling. Therefore it seems that the 
RWA can be applied for solving the JT problems. Most of the JT Hamiltonians are more complicated than the Dicke 
Hamiltonian. At present, a few of them (i.e. E^d, E^e ) has been analyzed in the framework of quasi- exactly solvable 
problem d, or isolated exact solvability [loL Ull Il3 Hsf , both provide finite number of exact eigenvalues and 

eigenfunctions in the closed form. In this paper we devise a novel method for solving JT Hamiltonians, as well as other 
quantum optical Hamiltonians in the framework of RWA. It will be shown that the eigenvalues and the associated 
eigenfunctions can be obtained in the closed form when the coupling constant is smaller than the natural frequency 
of the oscillator. 

The method described here includes a part of the motivation provided by the existence of the connection between 
JT Hamiltonians and Dicke Model. Here we concentrate our attention to the solution of the E ^ e JT Hamiltonian. 
Its solution has been treated previously by many authors 0, Q, 0, 0, 0, 0| . We develop a new approximation 
method which is based on the similarity transformation. The method introduced here is the same as the RWA which 
has been usually used to solve Dicke Hamiltonian. An interesting and somewhat simpler form of the JT Hamiltonian 
is obtained by RWA. 

Other purpose of this paper is to show that for some purely imaginary couplings the E®e JT Hamiltonian becomes 
non-Hermitian but its low-lying part of the spectrum is real. It will be shown that the non-Hermitian Hamiltonian is 
not PT-invariant [UpO, 21, 22, 23] , but it is pseudo-Hermitian [ilEil El EtL . 

In the following section, we shall demonstrate our procedure on the E ® e JT Hamiltonian. We present a trans- 
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formation procedure and we obtain approximate form of the E ® e JT Hamiltonian. We show that Hamiltonian can 
be transformed in the form of the Dicke type Hamiltonians. We also obtain exphcit expressions for the eigenstates 
and eigenvalues of the JT Hamiltonian. In section 3, we discuss the pseudo-Hermicity of the Hamiltonian. Finally we 
summarize our results. 



METHOD AND SUMMARY OF THE PREVIOUS RESULTS 



The well-known form of the E®e JT Hamiltonian describing a two- level fermionic subsystem coupled to two boson 
modes has obtained by Reikjl2| is given by 



H = oj (a+fli -I- aja2 + l) + t^o^o + ^[(ai + (^2)'^+ + i^i + 02)0'-], 



(1) 



where uq is the level separation, lu is the frequency of the oscillator and k is the coupling strength. The Pauli matrices 
tTo,± are given by 
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The annihilation and creation operators, and a^, satisfy the usual commutation relations, 



[fti, aj] = 0, [ai, aj j - ^..y . 



(3) 

The Hamiltonian |^ can be solved in the framework of quasi-exactly solvable problems or by using numerical 
diagonalization method^- In order to obtain rotating wave approximated form of the E ® e Hamiltonian, we use 
similarity transformation by introducing the operator 



T = 



cr-flo — (T+02] 



(4) 



and imposing the condition [w ± wqI 3> k, which usually holds in the weak interaction, transformation of the Hamil- 
tonian (0), yields that 



H 



He ^ K LO {alai + a^02 + l) + t^oO'o + ^[(ai + a2)(T+ + {a^ + aj)a_] + 



{0^02 + 0,102) + 



2 

to / -4-2 2 I \ 

"2 2 ("2 + ^2 + 02) O-Q 



O 



cifo2 + flia^ 



CTO 



(5) 



Since ^^^^ <C 1 is assumed to be a small parameter, neglection of the last term confirms result; 



(6) 



It is analytically solvable due to the neglect of the counter-rotating terms, so called RWA. Now, we turn our attention 
to the solution of the Hamiltonian ((HJ. The rotation of the bosons given by the following operator 



provides the expressions 



U = exp y — {o]\a2 — a2 fli 



U{afoi + 0202)1/^^ — OiOi + a^a2 



(7) 



Under [/, the Hamiltonian becomes 

H K UJ (o^ai + 0^02 + 1) + LUoao + V2K[aia+ + o'^a^]. 



(9) 
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Ground state 


First excited state 






Eexact 


Erwa 


Eexact 


0.1 


0.90455 


0.90442 


1.85982 


1.82286 


0.2 


0.81678 


0.81595 


1.73508 


1.67515 


0.3 


0.73508 


0.73277 


1.62159 


1.54472 


0.4 


0.65835 


0.65371 


1.51676 


1.36373 


0.5 


0.58578 


0.57798 


1.41886 


1.31592 


0.6 


0.51676 


0.50498 


1.32667 


1.21248 


0.7 


0.45080 


0.43429 


1.23931 


1.11438 


0.8 


0.38754 


0.36557 


1.15609 


1.02070 


0.9 


0.32667 


0.29856 


1.07646 


0.93072 



TABLE I; Ground-state and first excited-state energies of the E ® e JT Hamiltonian. 

The resultant Hamiltonian can easily be solved, because the matrix of the Hamiltonian can be decomposed in infinite 
dependent 2x2 blocks on the subspaces {|t, ni) , 712) , | J,, + 1) |i, ^2)}, where ni and 712 are the number of photons. 
The eigenvalue problem can be written as 



where \^) is the two component eigenstate 



H\^)=E\^) 



ci \ni) I712) 
C2 Ini -f 1) \n2) 



(10) 

(11) 

(12a) 
(12b) 



where ciand C2 are normalization constant. Action of H on -ip yields the following expressions 
(^ci {to {ni + 712 + 1) + t^o) + C2\/2K\/ni + 1) \n2) = ciE \ni) \n2) 
(^C2 {lu {ni + n2 -f 2) - luq) + ciV^Ky/ni + ij \ni + 1) |ri2) = C2E \ni + 1) |n2) . 
Eliminating ci and C2 between (|12a|) and (|12b|l and solving the resultant equation for E, we obtain 

E= {j + l)oj±^^j8K^{n + l) + {Lj^2ojaf. (13) 

, 2j. The eigenstates can be easily written by using 



where j = rii + n2 total number of bosons and n — 0,1, 2, 
boson operators, acting on a vacuum state |0) ; 



Cia2 



|0),C2arX"+ |0) 



(14) 



We conclude that in weak coupling limit the oscillators does not coupled to each other and each of them oscillates 
with their own frequencies. We have proven that when the interaction between E ion and e-modes are weak then 
E ® £ JT Hamiltonian can be reduced to the JC model. Our formalism provides a solution of the problem which 
allows us to discuss the JT effects in the Dicke model. 

The accuracy of the approximate eigenvalues can be checked by means of the (quasi) exact solution of the E ® e 
JT Hamiltonian. The material parameters are chosen to be w = 1 and = 0. The results are tabulated in Table 1. 

The results of our study show that the eigenvalues and eigenstates of the E®£ JT Hamiltonian can be approximately 
described when the frequency lo of the oscillation larger than the interaction constant. 



NON-HERMITIAN INTERACTION 



It has been shown that for some purely imaginary couplings constant k, the low-lying part of the E ® e JT 
Hamiltonian is real, although the Hamiltonian is non-Hermitian. Let us consider the Hamiltonian with the 
imaginary coupling k = ij : 



: w [a^ai + aja2 + l) + ^o^o + iV2j[aicr+ + a+a-_]. 



(15) 
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This Hamiltonian is not Hermitian as, 

= LO (flj^ai + a^a2 + l) + loqcto — iV2^[aia^ + afa-] ^ h. (16) 

Under the parity transformation, the Pauh matrices become invariant but both the creation and annihilation operators 
change sign. The time reversal operator for this Hamiltonian is T = —lUyK where K is complex conjugation operator. 
The time reversal operator changes the sign of the Pauli matrices and boson operators. It is easy to see that the 
Hamiltonian H15|l is not PT-symmetric 

{PT)h{PT)-^ = Lu {a+ai + a+aa + l) - woctq + iV2'y[aia+ + a^a-] ^ h. (17) 

The Hamiltonian is not PT-symmetric but it gives real spectrum. Mustafazadeh [j^. Esl |23| has shown that the 
reality of the spectrum of non-Hermitian Hamiltonian is due to pseudo-Hermicity properties of the Hamiltonian. A 
Hamiltonian is called 77-pseudo-Hermitian if it satisfies the following relation 

ry/iry-i = h\ (18) 

where is a linear Hermitian operator. The Hamiltonian h and its adjoint /i^ can be related to each others by the 
operator ctq and using the relation a^a^a^^ — —a± : 

G^ha^^ = h^. (19) 

Then the Hamiltonian 1)15(1 is cTQ-pseudo-Hermitian. Our Hamiltonian is also pseudo-Hermitian with respect to the 
parity operator. As it is shown 24] that if a Hamiltonian is pseudo-Hermitian under two different operators, r]i, 772 
then the system is symmetric under the transformation generated by 771 ?/2^^. Therefore our Hamiltonian is invariant 
under the symmetry generated by the combined operator, Pctq : 

[H, Pao] = 0. (20) 



CONCLUSION 



The aim of the this paper was to illustrate how the E ® e JT Hamiltonians can be solved by developing a trans- 
formation procedure. It has been found an approximate form of the E ® e JT Hamiltonian in the framework of the 
RWA. The resultant Hamiltonian can be solved analytically and its eigenvalues can be obtained in the closed form. 
We have shown that in the weak coupling limit the JT models may be recognized as the Dicke model. We have shown 
that when the coupling constant is imaginary the Hamiltonian is non-Hermitian but Pdo-symmetric. We also hope 
to extend the method to the other JT and quantum optical systems. 
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